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Abstract 

We study mesoscopic pairing in the one dimensional repulsive Hubbard model and its interplay 
with the BCS model in the canonical ensemble. The key tool is comparing the Bethe ansatz 
equations of the two models in the limit of small Coulomb repulsion. For the ordinary Hubbard 
interaction the BCS Bethe equations with infinite pairing coupling are recovered; a finite pairing 
is obtained by considering a further density-dependent phase-correlation in the hopping amplitude 
of the Hubbard model. We find that spin degrees of freedom in the Hubbard ground state are 
arranged in a state of the BCS type, where the Cooper-pairs form an non-condensed liquid on 
a "lattice" of single particle energies provided by the Hubbard charge degrees of freedom; the 
condensation in the BCS ground state corresponds to Hubbard excitations constituted by a sea of 
spin singlets. 

PACS numbers: 74.78.Na,04.20.Jb 
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I. INTRODUCTION 



Intensive analysis has been devoted to trace pairing phenomena in strongly interacting 
electronic model o^i^i^ . Most of the motivations come from studies in high-Tc superconduc- 
tivity for which it is highly desirable that pairing mechanisms are automatically encoded 
in the interaction. The Hubbard model constitutes the prototype of the class of models of 
interest in the subject^. It describes a lattice gas of itinerant electrons experiencing an on- 
site Coulomb interaction. So far, the evidence of pair-pair correlation is ambiguous both for 
repulsive and attractive Hubbard interaction^'^. Here we focus on the Hubbard (Id) chains, 
with repulsive interaction^. Not only does a Id analysis provide a unique theoretical lab 
to predict the behavior of compounds in higher dimension, but current technology actually 
admits quasi- Id superconductors to be fabricated; finally, there is evidence that the stripe 
ordered domains gives effectively a Id character to the phenomenology of even many su- 
perconducting 3(i-compounds^. However, how the superconductivity does arise from stripes 
still constitutes an open problem. Assuming that the electronic properties of the stripe- 
domains can be captured by repulsive Hubbard-like models, effective pair-formation can 
be studied by means of the pair binding energy Epf,, detecting an energy gain whenever 
an electron-electron attraction is established into the systen>^»i2iii. In particular, Eph was 
studied for the ground state (GS) of the finite Id Hubbard model and it turned out that 
the repulsive local interaction can be over-screened, leaving an effective attraction among 
the electrons^^. This phenomenon was called mesoscopic pairing. An important motivation 
for us is to single out a pairing of the BCS type in the mesoscopic pairing of the Hubbard 
model. The attempt to trace BCS pairing within the Hubbard type models is not new, and 
it has been well studied in the physical community, mainly for attractive interaction. A 
common feature of such studies^- (with the exception of the variational study in Ref. Il2l ) 
is that the analysis was done in the grand-canonical ensemble, involving mean field BCS 
(that is exact in thermodynamic limilji^). The scenario changes in the canonical ensemble 
where the BCS gap vanishes identically, the "superconducting phase" is dominated by quan- 
tum fluctuationsi^ and the mean field BCS theory is inadequatoi^. On the other hand, the 
canonical analysis is a natural tool to study the mesoscopic BCS pairing in the Hubbard 
model. This choice in the context of the stripe physics is supported by various experimental 
evidences. Stripe ordered domains in superconducting nano-powders were detected in Ref. 
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16l : fragmented nano-stripes (10 — 20 nanometers long) were clearly noticed also in bulk 
high-T superconductors^-. We use the Bethe ansatz solution of the BCS model to compare 
the two pictures in the canonical ensemble. The exact solution of the BCS model was found 
by RichardsoniSi in 1963, but went unnoticed to the condensed matter community. It was 
discovered only recently to study the superconductivity at nanoscal o^^i^^ . 

The paper is organized in the following manner. In the next section we introduce the Hub- 
bard model and its Bethe Ansatz due to Lieb and Wu. In section lTTTl a formal correspondence 
is established for [/ — > 0+, where the spin degrees of freedom (the spin "quasi-momenta" ) 
in the Hubbard GS can be seen as Cooper-pair quasi-energies corresponding to certain exci- 
tations of the BCS model. A finite pairing coupling can be achieved by considering certain 
gauge-field correlated- hopping in the Hubbard model, as shown in the section HVI The corre- 
spondence of the BCS GS to a "sea" of spin singlets on top of the Hubbard GS, is discussed 
in section The role of the BCS pairing in the mesoscopic pair-binding is discerned in the 
section EH Finally, section I VI II is devoted to our conclusions, and the discussion of some 
physical implications. 

II. THE HUBBARD MODEL 

We consider a Hubbard chain of length L, with N = N^ + Ni particles M = Ni spin down, 
and periodic boundary conditions. Formulated in second quantization, the Hamiltonian is 

L L 

H = -t Xl(4+i,<xCj,a + h.c.) + UJ2 ^jA^jd ' (1) 

where {cj>,c[^,} = 6a,a'Sj^i, {cj^a,ciy} = 0, and rii^^ ■= cla'^i^^ ^ = {T' i} is the electronic 
spin; U and t (we set t = 1) are the Coulomb interaction and the hopping amplitude 
respectively. The many-body wave function 

^ Yl I]^(-^'i'---'-?'^l^)sW)---4<x.(iv)|0)' (2) 



i<ii<---<ijv<i Sn 
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G(fl)= n "T (6) 



brings Eq. (|T]) to the first quantized spectral problem Tiip = Sip, that can be solved by the 
Bethe ansata^^ 

V'(ji,---,j7vk) = ^ A^{Q)exp <i^kQ^i)ji\ , (3) 

Q&Sn I 1=1 J 

^^(Q) = Yl G{R)l[F{XR(^^j,y^\Q), (4) 
2i[/ ^ A — sin — it/ 

F{X,yQ) = —r—^U-. . ' 5 

X — sm ky + tU A — sm ki + lU 
A« — Afl — 2iU 

l<a<l3<M 

In these equations vr labels a permutation of the indices in the configuration space, whereas 
Q and R label permutations in the {k}, {A} values in the symmetric group Sjy and Sm 
respectively. 

The exact energy and momentum of the model^*^^ 

N N 

E = -2Y,cosk, P = Y.^, (7) 

are given in terms of charge-momenta kj, solutions of the Lieb-Wu Bethe Equations (BE)^ 

^ sin i-. _ A + 



n sm fej - Aa + — ^ T-r Ag - A/3 - — 
dr, t X J- J- X X. -U ' ^ 

sm fcj - Aa- — p=i -^a - AfS + — 



where j G {1, 2, . . . A^} and a G {1, 2, . . . M}. 

The equations ©, (jH)), and constitute the exact solution of the Hubbard model both 
for repulsive and attractive interactions. We focus on the limit U ^ 0^. The result of such 
operation on Eqs. ((Zj), (jHl) and (jHI) depends on the actual configurations of the charge and 
spin rapidities k, X. For numerics, systems up to = 100 electrons have been considered. 

III. HUBBARD GROUND STATE FOR VANISHING INTERACTION 

Here we consider the half-filled chain, and = 2M, even. For U > the GS 
is characterized by real k and A, where any A^, is trapped between two sin ka- Systems 



with N = 4n electrons are characterized by degenerate sink = (see the central cluster 
of Fig. ^ (a)). At small U the two sinfca's are arranged symmetrically around as 
I A,- — sin ki I ~ CjU^/'^ + DjU where Cj and Dj are numerical coefficients; from the numerical 
analysis we find Dj = 1/4. We further assume the quasi- momenta having the following 
behavior kj ~ kj^^ + f/^/^/cj^^^'* + Uk^p; this assumption will be confirmed by our numerical 
studies. The limit [/ 0+ of Eqs. © and Q readsiSiSi 

e'^f"- = 1 , (10) 

^'^-hY.'^^^ (11) 

^ ^ sm k)j - Aq 
^ A/3 - Aa ^ sin kj -\a~^' ^^^^ 



and k^"^^ L = l/2Cj. Yl' Eq. (jllj) indicates that the term sinfco, is omitted. Notice that 
U has to be small in order that both lU^^^k^'^^ L\ <^ 1 and {Uk^ L\ -C 1 are satisfied and 
that degenerate sinfcj may occur in Eq. ()12|) (see Fig. Q). 

The Bethe ansatz wavefunction in the limit of vanishing U is^ 

(^ = n - sin[A:Q(,„)] ^^^^ 

The relevance of this limit resides in the following observation. The Eqs. ()12|) are known 
as the Gaudin equations. They play an important role in the exact solution of the BCS 
model at infinite pairing coupling g. The BCS Hamiltonian for Np time-reversed pairs in Q 
single particle levels ej with pair degeneracy dj reads^^ 

n n 

^ '"im^ifh^j m' rri' 



m j=l i)j=^ m,in' 

n 



-.-Errf + k^ (14) 



where c^,^ creates an electron in the state \j;Sj,m), m G {1 . . .dj}; \j; Sj,Th) is the time 
reversed state of \j]Sj,m). The operators S~ := J2m'^jrnCjm, Sj' := (5*^")^ 5"! : = 
iYlm=-dj ''^jm ~ dj)/2 form a dj/2 dimensional representation of the su{2) algebra. The 
second line of Eq. (fT^ estabhshes the relation between Hbcs and the twisted six-vertex 
models22, through the Gaudin model^^, ultimately proving the integrability of the BCS model. 
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FIG. 1: The spin (Aq, black dashed hnes) and charge {s'mkj, red sohd hnes) rapidities at U/t <C 1, 
for a system with 8 (a) and 10 (b) electrons. A correlated hopping ^ removes the degeneration 
of the cluster at A; = for = 4n (c). In the upper cluster, though exhibits a slight dependence 
on U, the functional behavior of the difference |Aa — sin/ca| remains the same (inset). 



Its exact energy and eigenstatesi^ are 



BCS 



a=l i=l 



1 



1^)- E n 



1 



R&Sa 



R(a) 



Np n 

HE 

Q = l j = l 



s] 

' -10) 



(15) 
(16) 
(17) 



where the second sum in Eq. fll5|l runs over the set of A^";, blocked levels, i.e. occupied by a 
single electron and not available to pair scattering. Different states in Eq. (fTTj) correspond 
to different distributions of the Cq, at = on the ^j-configuration; in the BCS GS all e/s 
are complex conjugates. The relation between Eqs. (fT ^ . (fT7j) and Eqs. (fT^. (fT^ is evident 
identifying 2ej and Cq, (the pairing quasi-energies) with sin kj and Aq, respectively, fl and Np 



playing the role of and M. Eqs. fjl2|) are obtained from Eq. p6p in the limit g ^ oo. 
Important for our purposes are the excitations of the BCS model characterized by real {cq,}; 
it was noticed that in the canonical ensemble these can be obtained as the states in Eq. (|17p 
corresponding to {cq,}, finite solutions of the Gaudin equations^^. Therefore we conclude 
that the spin-rapidities (all of them are real and finite) in the Hubbard GS are arranged 
along the configuration of the BCS excited state displayed in the Fig. El 

IV. FINITE PAIRING 

Now we show that a finite value of g can be obtained by adding a further interaction 
in the original Hubbard model, in form of correlated-hopping. The correlated hopping we 
consider corresponds to 



in Eq. (^. The Hubbard model with the correlated hopping Eq. p8|) are of the Shastry- 
Schulz typ©25.. They were solved exactly for A, a obeying certain restrictions^^. It was 
demonstrated^fi that the effect of the correlated hopping is to twist the boundary conditions. 



^Ti(f^) = Et=i«i.'"(^)' ^Tt(^) = EjVm-i,m^i,'n(o') + ^m,m-i(fT) + A^_i,^+i ((x) . The BE 



of the Schulz-Shastry models are obtained from the Lieb-Wu BE, by multiplying with the 
factors e~**T and e*^*T-*i) the r.h.s. of the equations (jH)) and respectively. To obtain 
a finite g, it is assumed also that $0- ^'"a^ + f/^a"*. For —it < kj < n we find that the 
conditions = — ^^^^ = and |f/<l>||''| <^ 1 must be satisfied, in order to obtain the 
BE to zero order in U. The Eqs. fllOj) and then modify into 




(18) 



The boundary phases are $o- := (l>{cr)+9i[{o')N^a-+(^'[i{o'){Na- — I), with (f){a) = J2j=i 




(0) 



e 




and Eqs.dni) into 




(19) 
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with g = 1/(2$!^''). Thus the BE of the correlated hopping Hubbard model lead, in the 
hmit [/ -> 0+, to Richardson BE, Eq. In the same limit the Eq. ^ is 

E = — 2^ + const. 



a=l 

with the constant being UM{N - M +1)/L - 2^J^^coskf - (2f/<l>}^V^) Ejli sin fcf ^ 
That is: for small U, the energy of the Hubbard model coincides (up to constants) with the 
energy of the BCS model , see Eq. lfT^ . 

The excitation of the BCS model corresponding to the Hubbard GS is achieved by filling 
the levels e/s with a Cooper pair, leaving one empty level between two filled ones. In 
this way, the empty levels prevent the pairing parameters to form complex conjugate pairs. 
Such a BCS excitation is constituted by an non-condensed liquid of Cooper pairs (see the 
discussion in the next paragraph). The structure of these excitations is different for = 
4n + 2 and = 4n. In the first case the ej are all non degenerate: dj = 1, Vj; for = 4n, 
instead, the doubly degenerate Fermi energy correspond to a doubly degenerate Ej, as one 
can see in the Figs. QandEl). 

V. THE BCS GROUND STATE 

The trait of the condensation of the Cooper pairs in the BCS GS is the emergence 
of complex conjugate Cq, solution of Eq. (fTHjl . In fact by resorting certain electrostatic 
analogy2L2£ it can be proved that the thermodynamic limit of the Richardson's BE leads to 
the gap equation: 2G p{s)/ ^/{e — a){e — b) = 1 where p{e) is the density of single particle 
levels; the parameters a = Eq + iA and b = Eq — iA are the the end points of the arc where 



the solutions of Eqs. lfT^ are disposed (see FigOl see also Refs. |22||29|). Thus: the BCS gap 
(in the thermodynamic limit) is directly related to the maximum of the imaginary part of 
the CcT^iSSiSS. For large g all the Cq, are complex conjugated pairs (with a single exception if 
Np is odd); their real part is far below the lowest Ej (see e.g. Ref. 

To find the Hubbard eigenstate corresponding to the BCS GS we insert into Eqs. (jH]), 
the ea{= Aq) configuration obtained by solving Eq. (fT^ with certain configuration of 
sinkj^^; this fixes the kj, given the $o- we have used for Eq. ()19p with ^ 1. We 

found that all the k/s are real and non degenerate for non degenerate sin A;j'''*'s. (degenerate 
sin k^^^ correspond to degenerate kj). Numerically we found that the mapping can be reliably 



FIG. 2: BCS and Gaudin configurations leading to finite solutions at large pairing, for N = An 
(a) and N = An + 2 (b) and n = 2. The plot shows the behavior of the pairing parameters while 
increasing black curves, whereas the single particle energies remain constant, red lines. On left 
sides of the plots, the initial configuration of the pairing parameters (black circles) on the energy 
levels 2ej (red crosses) is shown; on the right sides the final configuration is displayed. The latter 
are the charge and spin configurations corresponding to the ground state of the Hubbard model. 
The larger cross in plot (a) indicates the double degenerate level 2e = 0. 

done for U < 0.03. It then turns out that sin/cj = sin/cj°^ + USf and | sin A; — A| 7^ is 
always satisfied. These conditions imply that equation (jl9|) characterizes also the excited 
Hubbard state under consideration. Summarizing, real kj & complex conjugate Aq, pairs, 
solutions of the Richardson BE, satisfy the Lieb-Wu equations. We conclude that the BCS 
GS corresponds to a sea of spin singlets (A-2 strings) in the Hubbard chain. Such a state can 
be seen as a BCS condensation of spin rapidities. As an example of such result, the table 
shows the k, A solution of both the BCS and the Hubbard GS BE for U = 0.001, = 6, 
and M = 3. 



sin kj 


Aq 


-2.14679391, -1.09959633 


-14.8060834 


-0.052399128, 0.994798057 


-11.722088 - i 11.1576888 


2.04199559, 3.08919348 


-11.722088 + i 11.1576888 
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FIG. 3: (a) The real and (b) the imaginary part of the pairing parameters are shown as functions 
of the pairing strength g for a system with = 16 electrons at half-filling in the BCS ground state 
configuration. In the thermodynamic limit the pairing parameters, being all complex conjugate at 
sufficiently large pairing, are arranged along arcs (see Ref. l29l ). where the distance between the 



end-points is the BCS gap 2A. Here we have set g = 3.0 and N = 200, with the single particle 
levels uniformly distributed in the range [-1,1]. 

VI. MESOSCOPIC PAIRING 



Here, we have a look at the pair binding energy Epi^ = 2E{N + 1) — E{N) — E{N + 2). 
For the Hubbard GS, Epf, shows "super-even effects"^*ii. Effective pair attraction Epf, > 
occurs for N = 4n {n is integer), and it is related with a vanishing GS spin gap = 
E{N)s=i - E{N)s=o for U 0+; for N = 4n + 2, Epf, < 0, corresponding to A5 > 
at small U. In Fig.Q we present Epf, for the Hubbard and for Shastry-Schulz models. 
We observe that the correlated hopping, Eq. (fTHj) . washes out the super-even effects (see 
FigE)); for = 4n + 2 and U > Uc Ep > (see also Ref. |lo|); at small we found 
that Uc oc A^. For larger N, Uc decreases (at = 50, Uc ~ 2.5). The phenomenon will 
be analyzed elsewhere. The formation of A5 can be interpreted within the Hubbard-BCS 
correspondence: the state with = 1 is obtained from 5 = by breaking a BCS pair 
and thereby blocking single particle levels. In this scenario the super-even effect arises (see 
also FiglSJ because breaking a pair in the state with A^ = 4n is energetically favorable. 
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FIG. 4: The super-even effects for the pair binding energy Epb (a) and the spin gap (b). Thick 
and hght hnes are results for periodic and twisted boundaries respectively. = 10: black solid 
line and red dotted-dashed line; = 12: blue dotted line and orange dashed line. The boundary 



twist is set to = = 0.2tt, = O.lvr = 2^>J^\ 



creating two blocked degenerate single particle levels, without kinetic energy extra cost. In 
contrast, for N = An + 2 the blocking of levels occurs at A; > kp- Away from half filling 
this is sufficient to cause an increase of the total energy; at half filling the levels are blocked 
in such a way that all the levels in the band are reshuffled, since breaking a pair at the 
Fermi energy causes the occupancy of the state at k = 0. We study the phenomenon by 
considering the pair correlation function \E'j = \uj\, where uj = {c^j m'^] m'^j mCj m) detecting 
the fiuctuational superconductivity in the canonical ensemble (in the thermodynamic limit 



becomes the BCS gap)^. Here we exploit the exact results achieved in Ref. where 
\I'j = \uj\ is calculated from a generating functiorn^i. 

Exclusively for systems with pair attraction (A^ = An with = 0) a dip at /c = is 
observed (FigEl); this reflects the corresponding small energy gap required to occupy the 
k = mode in the S = ^ S = 1 process. 



VII. CONCLUSIONS 



We proved that the exact solution of the BCS model can be obtained from the Bethe 
ansatz solution of the Hubbard model at vanishing U. At a formal level we comment 
that the amplitudes At,{Q) diagonalizing the 'spin sector' of the Hubbard model are, in 
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FIG. 5: The energy configuration (— 2cos /c) of the charge rapidities in the S = 0,1 ground states 
for N = 4n (a) and = 4n + 2 (b) systems. The different symbols indicate the degeneracy of the 
levels: circle, cross and triangle stand for one, two and three-fold degenerate levels respectively. 
Note that the three-fold degenerate levels are more precisely "quasi-degenerate" ones: the energy 
difference cannot be resolved in the plot. At = 4n, passing from 5" = to S" = 1 state two 
degenerate single particle levels are created and no additional energy is required (as shown by the 
arrows). Instead, at = 4n + 2, at half filling the levels are reshuffled, a pair is broken at Fermi 
level and a A; = level is occupied (upper arrow). 




sin k. 

J 



FIG. 6: The pairing correlator = \uj\ is presented (in the thermodynamic limit would 
become the BCS gap). We notice the anomaly for N = 4n. 
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fact, eigenvectors of the transfer matrix of certain inhomogeneous six vertex modelsi^. The 
relation between the Hubbard chain and the BCS model arises because the BCS model is 
itself a quasi-classical descendent of the inhomogeneous six vertex model^^. The Coulomb 
repulsion U /t plays the role of the quasi-classical parameter; the charge degrees of freedom 
{sin A;} of the Hubbard model play the role of the inhomogeneities for the vertex models. 
In the BCS picture, they provide the "lattice" of single particle energies Sj where the spin 
degrees of freedom, the spin "quasi-momenta" Aq,, can condense. The Hubbard ground state 
corresponds to a certain excited BCS state (see Fig|21). We noticed that also the Hubbard 
low lying spin excitations corresponds to states of the BCS type. In particular we have 
demonstrated that the BCS ground state of the spin rapidities reflects the formation of a 
sea of spin singlets in the Hubbard chain. In this scenario the super-even effective attraction 
in the Hubbard and in the Shastry-Schulz ground states can be understood in terms of pair- 
breaking excitations in the BCS model. We have demonstrated that BCS-correlators serve 
to study such parity effects. 

The results obtained in this paper might be relevant for the analysis of the stripe order 
in the high-Tc compounds. 
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